Prospects of medium tomography using hard processes inside a soft medium 
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Hard processes leading to high transverse momentum hadron production are calculable in per- 
turbative Quantum Chromodynamics (pQCD) for proton-proton collisions. In heavy-ion collisions, 
such processes occur as well, and due to a separation of scales are calculable on the level of the 
production of high pr partons. These subsequently interact with the medium which originates 
from semi-hard and soft processes in the same collision before they hadronize. The presence of the 
medium thus modifies the momentum spectrum of outgoing partons and thus hard hadrons. The 
idea of jet tomography is to study this modification in order to learn about the medium density 
distribution. While this is a compelling idea, its practical application encounters some difficulties. 
In this paper, the capability of various observables to provide tomographic information about the 
medium is discussed and compared. 

PACS numbers: 25.75.-q,25.75.Gz 



I. INTRODUCTION 

The expression 'jet tomography' often used to describe 
the analysis of hard pQCD processes taking place inside 
the soft matter created in an ultrarelativistic heavy-ion 
collision. In particular the focus is on the nuclear sup- 
pression of hard hadrons in A-A collisions compared with 
the scaled expectation from p-p collisions due to loss of 
energy from the hard parton by interactions with the soft 
medium (see e.g. [H 0, 0, 0, ) • As the use of the word 
'tomography' suggests, the ultimate hope here is to make 
a measurement of the density profile of the soft medium. 
The main underlying idea of jet tomography is to start 
with a hard process which is believed to be well under- 
stood and calculable in vacuum. Such a process may 
be the production of high transverse momentum (pt) 
hadrons h, which in p-p collisions can be obtained in 
leading order perturbative QCD (LO pQCD) from the 
(schematical) convolution 



da NN^ h+ x = f i/N ( Xl ,Q 2 ) ® f j/N (x 2 , Q 2 ) 

fijk (1) 

where fi/N{ x iQ 2 ) are the non-perturbative distribution 
functions of a parton i to carry a momentum fraction x of 
the nucleon N when the distribution is resolved at a scale 
Q 2 (see e.g. [1,0])j ^ij^f+ k are tne elementary cross sec- 
tions for partonic scattering ij — > fk calculable in pQCD 
and Dj^ h {z, fiV) is the non-perturbative fragmentation 
function describing the hadronization of a parton / into 
a hadron h which carries the momentum fraction z of the 
original parton when the hadronization scale is /i 2 (see 
8, 9]). 
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If this process is placed into a medium, in principle all of 
these ingredients may be modified. It has to be remem- 
bered that the medium is created at the same time as the 
hard parton pair by semi-hard and soft inelastic scatter- 
ing processes between the nucleons of the colliding ions. 
However, if the momentum transfer in the hard partonic 
scattering is sufficiently high, then the separation in scale 
between hard scattering vertex and typical soft medium 
scale can be used to argue that &ij^f + k should not be 
influenced by the presence of a medium as the spatial 
resolution scale of the hard process is too small to detect 
medium constituents. 

On the other hand, it is well known that the parton distri- 
bution inside a nucleon being part of a nucleus is distinct 
from fi/N(x,Q 2 ) discussed above [HI EH 0] • However, 
the change of the distribution can be inferred from e-A 
or p-A scattering experiments and can hence be assumed 
to be sufficiently known in heavy- ion collisions. 
The medium is then first relevant for the parton emerg- 
ing from the hard vertex. In principle, one can as- 
sume that the presence of the medium modifies the 
hadronization process of the parton and introduce a 
medium-modification of the fragmentation function (see 
e.g. [l3|, It is experimentally well known that 

hadronization processes other than fragmentation, pos- 
sibly coalescence [l5| or parton recombination [n| domi- 
nate hadron production below 6 GeV transverse momen- 
tum. However, when the hadron is created at sufficiently 
high pt it can be assumed that hadronization takes place 
outside the medium. A rough estimate for the formation 
distance of a hadron is given by Pr/^had which for a 10 
GeV pion evaluates to 0(100) fm, i.e. the assumption 
that the hadronization is not influenced by the presence 
of a medium is rather safe. 

Under these conditions, the medium effect on the pro- 
cess can be cast into the form of a probability density 
P(AE) of the parton experiencing the energy loss AE 
while traversing the medium. This distribution can be 
dependent on the parton color charge (i.e. distinct for 
quarks and gluons) and can in principle also depend on 
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the initial parton energy. 

Thus, in the presence of a medium, for sufficiently hard 
processes the production of high pr hadrons can be writ- 
ten as 



/ 

(2) 

where 



dv± A r f+X = J2 fi/A&uQ 2 ) ® fj/A(x 2 , Q 2 ) ® 

ijfe 

(3) 

In these expressions, fi/A( x i> Q 2 ) now denote the nuclear 
parton distribution functions [10, EH and Pf(AE,E) is 
the probability of energy loss to the medium. In this kine- 
matic limit, the tomographic information on the medium 
is therefore contained in Pf(AE,E). A measurement of 
high pt hadron production can provide information to 
constrain this quantity. However, in order to obtain to- 
mographic information of the medium, it has to be real- 
ized that Pf(AE,E) depends on both the precise nature 
of the interaction of parton and medium and on the den- 
sity distribution of the medium. Thus, in order to obtain 
medium density information, the energy transfer from 
parton to medium must be assumed to be known. 
In the following, we will elaborate on the concepts out- 
lined above. First, we will demonstrate that the energy 
loss probability is not well constrained by the suppression 
seen in single hadron spectra. We will then turn to more 
differential measurements and show what additional in- 
formation can be obtained to better constrain the energy 
loss probability distribution and what can be concluded 
with regard to the density distribution if a specific model 
for energy loss is assumed. 



II. THE NUCLEAR SUPPRESSION FACTOR 

The most commonly studied observable in the context of 
parton energy loss is the nuclear suppression factor 



Raa(pt,v) 



d 2 N AA /dp T dy 
T AA (0)d 2 a NN /d PT dy- 



(4) 



This quantity is unity if an A-A collision can be viewed 
as a superposition of independent p-p collisions (i.e. no 
medium is created). Any deviation from unity indicates 
a non-trivial effect which has to be attributed to either 
energy loss or the nuclear parton distributions. It has 
been shown that the effect of the latter is small [l7[ ■ On 



the other hand, experimentally Raa for identified pions 
as a function of Pt seems to be flat out to 20 GeV and 
has a numerical value of about 0.2, i.e. the suppression 
of high pt hadrons is substantial 



Since Raa does not contain any spatial information, 
the production vertices of hard partons and their path 
through the medium have to be averaged out. Hard ver- 
tices (xo,yo) are distributed according to a probability 
density 



T A (ro + b/2)T A (ro - b/2) 
T AA (h) 



(5) 



where b is the impact parameter. The thickness func- 
tion is given by the nuclear density pa(y,z) as Ta(t) = 
J dzpA^J, z). Hence, if we can compute the energy loss 
probability distribution Pf{AE 1 E) pat h for a given path 
through the medium, the quantity to be folded with the 
pQCD parton spectrum and the fragmentation function 
is the average over all possible points of origin and paths, 
i.e. 



(P f (AE,E)) TA 

o2tt poo pOO 



1 

2^ 



• oo poo 

dx dy P(x ,yo)Pf(AE,E) path 

-OO J —OO 



(6) 



Any information about the medium density is now solely 
contained in Pf(AE, E) pat h- 



III. SENSITIVITY OF Raa 

Before we turn to a specific model of the interaction of 
the hard parton with the medium, let us discuss how well 
the functional form of (Pf(AE, E))t aa is constrained by 
current data. This can be done by inserting a trial ansatz 
for (P f {AE,E)) TAA into Eq. © [H]. 
In general, for a parton with given energy E, the effect 
of (Pf(AE, E))t aa may lead to the following three possi- 
bilities: 1) the parton escapes without energy loss (trans- 
mission) 2) the parton experiences so much energy loss 
that it is shifted to a thermal range, i.e. it can no longer 
be considered hard but becomes part of the bulk medium 
(absorption) and 3) the parton momentum is shifted, but 
it still remains sufficiently hard after energy loss (shift). 
The definition of what AE leads to shift and what is 
absorption naturally depends on the initial parton en- 
ergy E. This is an implicit dependence of the energy loss 
probability which is relevant even if Pf (AE) pat h is itself 
independent of E. In addition, there may also be an ex- 
plicit dependence of the functional form of the energy 
loss probability on the initial parton energy. 
Several trial energy loss probabilities are shown in Fig.[TJ 
left panel. The different scenarios are described in greater 
detail in [l9j . In each, one parameter has been tuned to 
result in the best possible description of Raa- In ad- 
dition, we also test a constant fractional energy loss in 
which every parton loses 30% of its initial energy 
The resulting Raa after evaluating Eqs. (|2l4p is shown 
in Fig. [TJ right panel. Considering that all scenarios 
must underpredict the data below 5-6 GeV momentum 
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FIG. 1: Left panel: Trial energy loss probability distributions {P(AE, E))t aa and the result of a full hydrodynamical simulation. 
Right: Nuclear suppression factor Raa for different toy models and the hydrodynamical simulation compared with PHENIX 
data [H|. 



as the model does not include hadronization mechanisms 
other than fragmentation which is known not to describe 
the data in this region, it is apparent that present data 
cannot distinguish well between the different scenarios. 
Given an assumption about the functional forum of the 
energy loss probability, the data can only constrain a 
scale. An exception to this is the fractional energy loss 
which seems to display the wrong qualitative trend. One 
may conclude that the explicit energy dependence of the 
energy loss probability distribution is therefore weak. 

In detail, there are of course difference between the sce- 
narios, reflecting the presence or absence of an explicit 
scale for the average energy loss and hence the implicit 
energy dependence of (P/(A£ I , E))t aa - The typical en- 
ergy loss shows the strongest variation as a function of pr 
as the definition of shift and absorption changes across 
different px whereas the geometrical energy loss in which 
partons are absorbed based on their region of origin is 
scale-blind and shows no variation. 



7 dN 
PT dp h T 



dzD™J h (z,Q') I df(P q (f,pl))6[zf 
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(7) 



typical energy loss 

smoothed geometrical suppression 

semi-opaque medium 

hydrodynamics 

geometrical suppression 



i 

i 

_i 

10 



P T [GeV] 



IV. 7-HADRON CORRELATIONS 



It is difficult to exploit the scale sensitivity as a tool to 
distinguish between different (Pf(AE, E))t aa from Raa 
as Eq. ^ contains a convolution over all possible initial 
parton energies. This averaging clearly removes sensi- 
tivity. However, by using 7-hadron correlation measure- 
ments it is possible to create a monochromatic source of 
quarks propagating through the medium [l9|, [2(| [2l[ . 

The momentum distribution of hadrons observed back to 
back with a 7 (neglecting intrinsic kr smearing) can be 
found from 



FIG. 2: Momentum spectrum of hard hadrons correlated 
back-to-back with a photon trigger, normalized to the vac- 
uum expectation. 

This expression is evaluated with the trial probability 
distributions shown in Fig. [1] and the result is shown in 
Fig. [2] As clearly seen, in such a measurement the dif- 
ferent functional forms of the energy loss distributions 
can be more easily distinguished. While this technique 
allows in principle a model-independent measurement of 
(P q (AE, E))t aa , it is experimentally challenging and the 
result still does not disentangle the physics of energy loss 
with the geometry of density distribution. Let us there- 
fore turn to a specific model of parton interaction with 



the medium and discuss what can be learned about the 
medium under this assumption. 



V. RADIATIVE ENERGY LOSS 

Our calculation follows the BDMPS formalism for radia- 
tive energy loss [22| using quenching weights as intro- 
duced by Salgado and Wiedemann [23]. Detailed infor- 
mation about the implementation is given in [24| , here 
we just present an overview. 

If we call the angle between outgoing parton and the re- 
action plane cj), the path of a given parton through the 
medium £(r) is specified by (i"o,</>) and we compute the 
energy loss probability Pf(AE) pat h for this path by eval- 
uating the line integrals 

/>oo />oo 

w c (r o ,0=/ dtiq(0 and (qL){r , </>)=/ 
Jo JO 

(8) 

along the path where we assume the relation 



0.75 



< 



0.5 



0.25 




A A v 

STAR R cp (pions, x 0.75) 
STAR R cp (p+pbar, x 0.75) 
R AA from AKK (pions) 
R,, from AKK (p+pbar) 



feed 



10 

p T [GeV/c] 



15 



20 



FIG. 3: Raa for pions as measured by PHENIX [1SI | . scaled 
Rcp for pions measured by STAR [2jj (in good agreement 
with the PHENIX Raa), scaled Rcp for protons measured by 
STAR [2S| and calculated Raa for pions and protons, dashed 
lines indicate the regime where recombination is expected to 
become an important effect fit} ]. 



?(0 = K -2- e 3/4 (£)(coshp - sinhpcosa) (9) 



VI. Raa FOR IDENTIFIED BARYONS 



between the local transport coefficient (specifying 
the quenching power of the medium) and energy density 
e and the local flow velocity p where a is the angle be- 
tween parton propagation and flow direction. Here, oj c 
is the characteristic gluon frequency, setting the scale of 
the energy loss probability distribution and (qL) is a mea- 
sure of the path-length, weighted by the local quenching 
power. We view the parameter K as a tool to account 
for the uncertainty in the selection of a s and possible 
non-perturbative effects increasing the quenching power 
of the medium (see discussion in [241]). 
Using the numerical results of (23[, we obtain Pf(AE) pat h 
for uj c and R — 2oj^/ (qL) as a function of jet produc- 
tion vertex and the angle <f> (here R is a dimensionless 
quantity needed as input for the energy loss probability 
distributions as defined in (23j|). The energy loss proba- 
bility Pf(A E) vat h is derived in the limit of infinite parton 
energy [23|. In order to account for the finite energy of 
the partons we truncate Pf(AE) at AE — E- ]ct and add 
S(A — -Ejet) Je deP(e). This procedure is known as non- 
reweighting [2^. We point out that the alternative con- 
cept of reweighting to our understanding systematically 
overestimates Pf(AE) for AE <C £j e t and should be dis- 
regarded. In fact for a dense medium increasing q and 
employing reweighting leads to an increased escape prob- 
ability whereas increasing q and non-reweighting leads to 
the expected decrease in escape probability, see also [26[ . 
The further strategy is then to decide on a particular 
model for the dynamical evolution of the medium and to 
adjust K such that Raa for central collisions is repro- 
duced (the values of K needed vary about 50% within 
the class of models describing bulk data, see [27J). We 
then look at more differential observables to study where 
information beyond the basic scale K can be found. 



Based on the observation that the older KKP fragmenta- 
tion function [8j underestimates proton production in p-p 
collisions measured by the STAR collaboration [29( while 
the AKK set @ describes the data better, it has been 
suggested that the production of protons and antipro- 
tons in hard processes is gluon-dominated. This gives 
rise to the expectation that proton suppression should 
be stronger than pion suppression, as the gluon is more 
strongly interacting with the medium due to its different 
color charge. 

Since the number of produced quarks and gluons at given 
Pt is fixed from the pQCD calculation, the only change 
of the model going from pion to proton production is the 
fragmentation function, i.e. there is no further parame- 
ter adjustment. Protonic Raa is therefore a non-trivial 
prediction of the model. 

In Fig. [3] we show the model calculation of Raa with 
the energy loss implemented as described above for a 3d 
hydrodynamical model for the soft matter evolution [3(| 
for both pion and proton production and compare with 
PHENIX and STAR data Since our framework, 

making use of thermodynamical quantities to describe 
energy loss induced by the soft medium is not suited to 
describe 60 — 80% peripheral collisions where a dilute, 
predominantly hadronic, medium is expected to exist, we 
calculate Raa and scale the experimental data. We show 
that the scaled pionic Rcp as obtained by STAR [28| is 
in perfect agreement with the pionic Raa as obtained by 
PHENIX [18|]. This gives us confidence that the com- 
parison of protonic Raa with scaled Rcp is not grossly 
unreasonable. As apparent from the results, the mea- 
sured similarity of the pionic and protonic suppression 
is expected from our model calculations, the difference 
between pion and proton Raa is less than the error on 
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FIG. 4: pr dependence of Raa in plane (solid) and out of 
plane (dashed) emission at different values of impact param- 
eter b. 

the pion Raa over a large kinematic range. We stress 
that we consider the calculation to be an upper limit 
for the expected difference as the use of the AKK set of 
fragmentation functions overpredicts proton production 
in p-p collisions 29]. If this baseline is adjusted by hand, 
the difference between pion and proton Raa is reduced 
even further by 50% (not shown). 

VII. Raa VS. REACTION PLANE 

One can also go to a more differential probe by determin- 
ing the reaction plane in non-central collisions and mea- 
suring the suppression pattern of hadrons as a function 
of the angle between hadron and reaction plane. Since 
the initial overlap zone of the colliding nuclei is almond- 
shaped, hadron emission can be in-plane (i.e. probing 
the short distance of the almond) or out-of-plane (along 
the long axis of the almond). 

In order to calculate this quantity, the angular averag- 
ing in Eq. ([6]) has to be removed, in this way, the Raa 
acquires an explicit dependence on the angle <fi [27| . 
In Fig. @] we show the resulting Raa for emission in plane 
and out of plane for three different impact parameters for 
the medium specified by the 3d hydrodynamics [13] ■ The 
observed trend is as expected — the overall quenching 
power decreases with increased impact parameter as the 
produced medium becomes less and less dense. Further- 
more, an asymmetry between in-plane and out-of-plane 
emission develops with increased impact parameter, re- 
flecting the increased spatial asymmetry of the initial 
density distribution. 

However, it is well known that the pressure of the 
initial state creates elliptic flow which maps the ini- 
tial anisotropy in position space into a corresponding 
anisotropy in momentum space which removes the spatial 
asymmetry over time. The hydrodynamical model for the 
soft medium describes this feature of bulk dynamics well 
[30j , Thus, one may ask if the fact that the asymmetry is 
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FIG. 5: Raa for in plane and out of plane emission as a func- 
tion of pr at b = 7.5 fm, assuming a static medium given by 
the hydrodynamical initial state, a static medium with read- 
justed quenching power given by the hydrodynamical initial 
state and the full hydrodynamical evolution. 

removed over time is actually reflected in the observable 
Raa- I n Fig. [5] we show a comparison of the full hy- 
drodynamical evolution and a scenario in which we keep 
the initial state fixed (in this case, one has to rescale the 
parameter K downwards to get the same description of 
Raa for central collisions as the medium does not expand 
and the density remains high) . If one makes the compari- 
son, one finds that elliptic flow tends to decrease the split 
between in-plane and out-of-plane emission as expected. 
In addition, there is a small distortion in shape. 
A precise observation of Raa vs. reaction plane at suffi- 
ciently high pt, taken together with the requirement that 
the hydrodynamics must describe the measured elliptic 
flow, could therefore constrain the amount of elliptic flow 
which is built up early on (as jet quenching is dominated 
by the partonic evolution phase) vs. the amount of ellip- 
tic flow built up late during the hadronic evolution. 

VIII. DIHADRON CORRELATIONS 

Instead of studying the momentum spectrum of hadrons 
back-to-back with a 7, one can also investigate the spec- 
trum of associated hadrons back to back with another 
hadron. However, such a process shows pronounced dif- 
ferences. 

Experimentally, a hard hadron in a given momentum 
range is used as the trigger and correlations are ob- 
tained with respect to this hadron. However, due to 
the hadronization process, it is impossible to conclude 
from the trigger momentum to the underlying parton mo- 
mentum, rather the momentum smearing induced by the 
fragmentation function has to be accepted. 
But second and more important, since the trigger side 
(near side) parton is, unlike a 7, strongly interacting, it 
undergoes energy loss dependent on the amount of mat- 
ter it traverses. Thus, the probability density of find- 
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ing vertices in the transverse plane which lead to a trig- 
gered event is biased towards the surface of the medium. 
This in turn implies that when calculating the distribu- 
tion of partons back-to-back with the trigger (away side), 
Pf(AE) pa th is not averaged using Eq. ([6]) but rather over 
a conditional probability distribution which becomes ev- 
ident only after one has evaluated the trigger condition. 
Since this conditional probability density is nontrivial, it 
contains additional information on the medium density 
distribution. 

We use a Monte-Carlo implementation to simulate the 
problem. First, we sample the distribution of partons 
emerging from a hard vertex determined by Eq. (|3|). This 
yields the parton type (quark or gluon) as well as the 
transverse momentum. We define randomly one of the 
partons as 'near side' and propagate it to the surface of 
the medium. Along the path, we determine tu c and (qL) 
by evaluating Eq. ([5]). The resulting values again serve 
as input for the probability distribution of energy loss 
P(AE) pat h as determined in (23j . 

We determine the actual energy loss of the near side 
parton by sampling P{AE) pat h- To find the energy of 
the leading hadron, we need the probability P/^/j(z,/x) 
to find a leading hadron with momentum fraction z at 
scale /x. Strictly speaking, this is not the fragmenta- 
tion function as the fragmentation function yields the full 
single hadron distribution, not only the leading hadron, 
but since the trigger condition enforces on average large 
values of z, the two are virtually identical and we use 
the (normalized) fragmentation function Dj°^ h (z, fi) as 
a model for Pf^^iz, pi). 

If the resulting hadronic Ph = zpf f» zEf fulfills the 
trigger condition we accept the event and proceed with 
the calculation of associated hadrons and the away side 
parton, otherwise we reject the event and continue the 
MC sampling by generating a new vertex. 
If an event fulfilling the trigger has been created, we de- 
termine the kx smearing being added to the away side 
parton momentum. We sample a Gaussian distribution 
chosen such that the widening of the away side cone with- 
out a medium is reproduced. Since this is a number of 
order 1 GeV whereas partons fulfilling trigger conditions 
have frequently in excess of 15 GeV we note that this is 
a small correction. 

We treat the far side parton exactly like the near side 
parton, i.e. we evaluate Eqs. © along the path and 
find the actual energy loss from P(AE) with u> c , (qL) as 
input. If the away side parton emerges with a finite en- 
ergy, we again use the normalized fragmentation function 
Df°fi h (zi [i) to determine the momentum of the leading 
away side hadron. If this momentum fulfills the imposed 
trigger condition for associated particle production, we 
count the event as 'punchthrough'. 

In addition, we allow for the possibility that the fragmen- 
tation of near and away side parton produces more than 
one hard hadron. We cannot simply subsume this in the 
fragmentation function as done for single hadron distri- 
butions as we are explicitly interested in the correlation 



strength between near side trigger hadron and other near 
side hadrons, thus we have to calculate subleading frag- 
mentation processes separately. The quantity we need is 
Pf^i(zi,Z2,fJ-), i.e. the conditional probability to find a 
hadron i from a parent parton / with momentum frac- 
tion Z2 given that we already produced a leading hadron h 
with momentum fraction z\. In this language, the whole 
jet arises from a tower of conditional probabilities for 
higher order fragmentation processes (for a different ap- 
proach see [32|). However, since we only probe the part 
of this tower resulting in hard hadrons, the treatment 
simplifies considerable. 

Moreover, since we are predominantly interested in the 
quenching properties of the medium and not in detailed 
modelling of hadron distributions inside the jet, we model 
the next-to-leading conditional fragmentation probabil- 
ity using the measured probability distribution Ai(zp) of 
associated hadron production in d-Au collisions [33 . l34j 
as a function of zt where zt is the fraction of the trig- 
ger hadron momentum carried by the associated hadron. 
We include a factor 9(Ei — trigger — AE — E assoc ) on the 
near side and 9(Ei — -E pun ch — AE — E assoc ) on the far 
side to make sure that energy is conserved. Note that 
associated production on the far side above the pt cut is 
only possible if a punchthrough occurs. We count these 
events as 'near side associate production' and 'away side 
associated production'. 

Thus, the yield per trigger for dihadron correlations on 
the near side is determined by the sum of all 'near side 
associate production' events divided by the number of 
events fulfilling the trigger, the yield per trigger on the 
away side is given by the sum of 'punchthrough' and 
'away side associated production' divided by the num- 
ber of events (fulfilling the near side trigger condition). 
These quantities can be directly compared to experiment. 
In [24], [35| we have investigated dihadron correlations six 
different scenarios for the spacetime evolution of bulk 
matter from a parametrized evolution model |36l| as well 
as a 2d hydrodynamical evolution model [371 ]. In the 
following, we will mainly discuss the 2d-hydrodynamical 
evolution ('hydro') where we assume that energy loss 
happens throughout the evolution and the same evolu- 
tion ('black core') in which we confine energy loss to 
the partonic phase (with a corresponding rescaling of 
the quenching power to very large K). The first sce- 
nario corresponds to an extended medium with moder- 
ate quenching properties, the second to a dense core with 
strong quenching power surrounded by a halo from which 
partons can escape essentially without significant energy 
loss. Detailed information on all scenarios can be found 
in H. 

We show the calculated yield per trigger compared with 
the experimental result in Fig. [6j Within errors, the near 
side yield per trigger is described by all the models well. 
There is no significant disagreement among the models. 
The model calculations appear considerably more differ- 
ent if we consider the away side yield. Here, results for 
the 4-6 GeV momentum bin differ by almost a factor two. 
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presentation. 



However, none of the model calculations describes the 
data in this bin. This is in fact not at all surprising as be- 
low 5 GeV the inclusive single hadron transverse momen- 
tum spectra are not dominated by pQCD fragmentation 
and energy losses but, rather, by hydrodynamics possibly 
supplemented with recombination [TBI, [li| type phenom- 
ena. For this reason, the ratio Raa at pr < 5 GeV cannot 
be expected to be described by pQCD fragmentation and 
energy losses, either. This is clearly unfortunate, as the 
model results are considerable closer to the experimental 
result in the 6+ momentum bin on the away side and 
hence our ability to discriminate between different mod- 
els is reduced. Since at this large transverse momenta 
the pQCD fragmentation + energy losses dominate the 
singe hadron spectrum we expect that the model is able 
to give a valid description of the relevant physics in this 
bin: Not only is Raa well described by the data, but also 
the contribution of recombination processes to the yield 
is expected to be small 16]. Thus, as it stands, only the 
black core scenario can be ruled out by the data. 

Nevertheless, while there are indications that the other 
scenarios show sizable differences in the momentum spec- 
trum of away side hadrons, the present data is not suf- 
ficient to make a distinction. There are in principle two 
ways to overcome this problem. At the price of intro- 
ducing additional model dependence, one might include 
recombination processes into the simulation. In this way, 
comparison to the 4-6 GeV momentum bin would be pos- 
sible. Alternatively, one can address the question if more 
leverage in px would improve the question and hence if 
improved experimental conditions will allow tomography. 
We have chosen to follow the latter path. Towards this 
end, we show in Fig. [7] a situation where the trigger mo- 
mentum is increased to 12 GeV < pr < 20 GeV and 
consequently more momentum bins in the pQCD region 
become accessible. Clearly, the difference between the 
different evolution scenarios is somewhat magnified, al- 
though still not large. 



IX. THE GEOMETRY OF DIHADRON 
CORRELATIONS 

As stated above, there is a systematical surface bias of 
the vertex distribution induced by the requirement to 
observe a trigger. For two scenarios, this distribution is 
shown in Fig. [5] 

It is evident that the strength of this surface bias depends 
on the distribution of the medium quenching power — 
the surface bias of the 'black core' scenario where the 
central region is almost impenetrable is much stronger 
than in the 'hydrodynamics' scenario. Given that the in- 
medium pathlength of the away side parton is sampled 
from this distribution (rather than from Eq. (|5J)) the dif- 
ferences in the momentum spectrum of correlated away 
side hadrons in the two scenarios can be understood. 
This becomes even more apparent when comparing the 
vertex distribution of the subclass of events in which hard 
hadrons were observed on both the near and away side. 
This is shown in Fig. [9j It is quite evident that the pres- 
ence of a strongly quenching core leads to tangential pro- 
cesses in which the core is never traversed whereas in a 
more extended medium with moderate quenching power 
as in the hydrodynamics scenario the whole volume is 
probed. 

While these figures seem compelling, it has to be kept in 
mind that they have been calculated under the assump- 
tion that a specific model for the interaction between 
hard parton and bulk medium is realized. In this sense, 
the calculation represents a test of particular combina- 
tions of energy loss model and medium evolution rather 
than a model-independent truth. From this point of view, 
we can note that the radiative energy loss model used 
here based on [22], [23| can describe the correlation data 
assuming a range of hydrodynamical and hydro-inspired 
evolution models as long as quenching is not confined to 
the partonic phase alone. At present, the data are in- 
sufficient to resolve differences in detail between model 
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FIG. 8: Probability density of finding a parton production vertex at (x, y) given a triggered event with 8 GeV < pr < 15 
GeV for different spacetime evolution scenarios (see text). In all cases the near side (triggered ) hadron propagates to the ~x 
direction, hence the y — (—y) symmetrization. Countours are at linear intervals. 



evolutions. 



X. CONCLUSIONS 

We have discussed jet tomography, i.e. the idea to obtain 
information about the medium produced in heavy-ion 
collisions by studying the modifications of hard processes 
taking place inside the medium. We have argued that the 
tomographical information is contained exclusively in the 
energy loss probability of a parton / along a specified 



path Pf (AE, E) pa th if the measured hadrons are hard 
enough such that fragmentation can be assumed to be the 
dominant mechanism of hadronization and hadronization 
can be assumed to happen outside the medium. 

Different hard observables probe different averages of this 
quantity. The nuclear suppression factor Raa for pi- 
ons probes the average over all angles and the initial 
nuclear overlap and all parton species. By looking at 
proton Raa, one may assume that predominantly gluon 
suppression is probed. In contrast, by using 7-hadron 
correlations one gets access to quark energy loss and in 
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FIG. 9: Probability density for finding a vertex at (a;, y) leading to a triggered event with 8 GeV < pr < 15 GeV and in 
addition an away side hadron with 4 GeV < pr < 6 GeV for different spacetime evolution scenarios (see text). In all cases the 
near side hadron propagates to the —x direction. Countours are at linear intervals. 



addition a handle on the initial energy of the parton. 
In Raa vs. reaction plane measurements, the averaging is 
still performed over the full initial overlap geometry, but 
not over all angles which provides a systematic handle 
on pathlength dependence. Furthermore, low pt elliptic 
flow is tied to the removal of the spatial asymmetry which 
is reflected in the angular dependence of Raa and we 
have demonstrated that there is some sensitivity to the 
timescale at which this asymmetry is removed. 
Finally, back-to-back hadron correlation introduce (due 
to the trigger bias) a different distribution of vertices 
which is moved towards the surface. This implies that the 
away side parton is systematically biased towards long 
pathlcngths. While difficult to control theoretically, the 
changed geometrical averaging implies that additional in- 
formation on medium properties can be deduced from the 
momentum spectrum of away side hadrons. 
It is quite clear that taken alone none of the mea- 
surements outlined here is able to provide a model- 
independent determination of the medium density. Fur- 
thermore, there is a general problem that in a situation in 
which spectra are steeply dropping as a function of px ob- 
servables will always be dominated by events in which no 
energy loss has taken place. For this reason, any effects 



which are sensitive to the medium density distribution 
are in general subleading. 

It seems that the best strategy is not to try to determine 
medium properties not by considering hard probes alone 
but to use hard probes to test models which have been 
used to describe soft bulk dynamics and remove possi- 
ble ambiguities this way. While a single hard observable 
does not carry too much information, the collection of 
all measurable quantities sums up to to rather stringent 
constraints. Clearly, with both a large kinematic lever- 
arm for hard probes at LHC and high-statistics data from 
RHIC, hard probes will reach the precision needed to do 
tomography in the true sense of the word soon. 
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